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Abstract 



Supersymmetric models with an inverted mass hierarchy (IMH: multi-TeV 
first and second generation matter scalars, and sub-TeV third generation and 
Higgs scalars) have been proposed to ameliorate phenomenological problems 
arising from flavor changing neutral currents (FCNCs) and CP violating pro- 
cesses, while satisfying conditions of naturalness. Models with an IMH already 
in place at the GUT scale have been shown to be constrained in that for many 
model parameter choices, the top squark squared mass is driven to negative 
values. We delineate regions of parameter space where viable models with a 
GUT scale IMH can be generated. We find that larger values of GUT scale 
first and second generation scalar masses act to suppress third generation 
scalars, leading to acceptable solutions if GUT scale gaugino masses are large 
enough. We show examples of viable models and comment on their character- 
istic features. For example, in these models the gluino mass is bounded from 
below, and effectively decouples, whilst third generation scalars remain at 
sub-TeV levels. While possibly fulfilling criteria of naturalness, these models 
present challenges for detection at future pp and e"'"e~ collider experiments. 

PACS numbers: 14.80.Ly, 13.85. Qk, 11.30.Pb 
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Supersymmetry offers an elegant solution to the problem of quadratically divergent scalar 
masses in the Standard Model (SM), provided supersymmetric matter exists at or near the 
weak scale A Minimal Supersymmetric Standard Model (MSSM) can be constructed, 
with 124 free parameters, most of which occur in the soft SUSY breaking (SSB) sector of 
the model [0] and reflect our ignorance about how SUSY is broken. Taking arbitrary weak 
scale SSB parameter choices generally leads to conflict with various low energy constraints 
associated with flavor changing neutral currents (FCNCs), and CP violating processes such 
as the electric dipole moments of the proton and neutron 0]. Of course, SUSY model builders 
have to explain the origin of SSB terms while at the same time satisfying constraints imposed 
by low energy processes. 

Three possibilities have emerged for building models consistent with low energy con- 
straints: 1. universality (degeneracy) of scalar masses 2. alignment of fermion and 
sfermion mass matrices and 3. decoupling, which basically involves setting sparticle 
masses to such high values that SUSY loop effects are suppressed relative to SM loops 0. 
Models with gauge mediation 0], anomaly mediation p or gaugino mediation of SUSY 
breaking naturally lead to universality of particles with the same gauge quantum numbers. 
Supersymmetric models with SUSY breaking communicated via gravity (supergravity mod- 
els) in general lead to non-universal scalar masses 110]. The minimal supergravity model 



(mSUGRA) adopts universality as an ad hoc assumption [|l|]. In this paper, we explore a 
class of models which potentially solve the SUSY flavor and CP problems via the decoupling 
solution. 



It is important to notice that "naturalness" arguments [|ri|], which generally require 
sub-TeV sparticle masses, most directly apply to third generation superpartners, owing to 
their large Yukawa couplings. In contrast, the constraints from flavor physics mentioned 
above apply (mainly) to scalar masses of just the first two generations. This observation 
has motivated the construction of a variety of models, collectively known as inverted mass 



hierarchy (IMH) models |T^, where the first and second generation squarks and sleptons 
have multi-TeV masses, while third generation scalars have sub-TeV masses. For models in 
which the IMH occurs at or near the GUT scale (GSIMH models), it has been emphasized 
p!3| that two loop contributions to renormalization group (RG) running can cause tachyonic 
third generation squark masses to occur, unless these masses are beyond ~ 1 TeV, which 
again pushes the model towards the "unnatural" . 

Recently, it has been pointed out that models with a weak scale IMH can be generated 
radiatively by starting with multi-TeV scalar masses for all scalars at Mqut 0|- For certain 
choices of GUT scale SSB boundary conditions, and assuming Yukawa coupling unification, 
the Higgs and third generation SSB masses then evolve rapidly towards zero, whilst first and 
second generation scalars remain heavy. However, requiring realistic third generation fermion 
masses and also a consistent radiative breakdown in electroweak symmetry (REWSB), only 
a rather small IMH can be generated |T^, which is not sufficient by itself to solve the SUSY 
flavor and CP problems. 

An alternative approach to a decoupling solution is to take very large values of scalar 
masses in models with intermediate to large values of tan (3. The "focus point" behavior of 
the Higgs SSB masses results in models with all matter scalar above a TeV, but with low 



values of and possibly low fine-tuning [IG]. However, even in these models, scalar masses 



are typically in the 1-3 TeV range, and are again not sufficient to solve the SUSY flavor and 
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CP problems without some degeneracy or alignment. 

In this paper, we examine models with a scalar IMH already in place at the GUT scale. 
We assume that the MSSM is a valid theory between Mqut and M^^ak, and that REWSB 



occurs. As noted in Ref. [|1^, it is crucial to work with two loop RGEs for this class of 



models. The two loop RGEs for the MSSM have been presented in [|1^], and have been 
implemented in ISAJET versions > 7.49 [0. 

We adopt the following parameter space for our studies: 

mo(l), mo(3), mi/2, ^o, tan/3, sign{fx), (1) 

where mo(l) is the common mass of all first generation matter scalars at Mqut, and mo (3) is 
the common third generation scalar mass. For simplicity, we adopt r?7,o(2) = mo(l), although 
the whole point is that mo(2) need not equal mo(l), so long as both are far above the TeV 
scale. The Higgs scalar SSB masses are set equal to mo (3), and take values of ~ 1 TeV. As 
usual, mi/2 and Aq are common GUT scale gaugino and trilinear masses, tan/3 = and /i 
is the superpotential Higgs mass term. 

For our renormalization group solution to the sparticle and Higgs mass spectrum, we use 
ISASUGRA (a part of the ISAJET package). Briefly, starting from weak scale values for 
the gauge and Yukawa couplings, ISASUGRA evolves the couplings up in energy until the 
GUT scale is determined, where gi = g2. At Mqut, the various SSB mass parameters are 
entered, and the set of 26 coupled RGEs for gauge and Yukawa couplings, and SSB masses, 
are evolved down to scale M^jeak, where the sparticle mass spectrum can be calculated. An 
iterative procedure is adopted so that the RG improved one loop effective potential can be 
calculated and minimized at an optimized scale choice Q = y/rn^^ m^^ , where REWSB is 
required. 



The form of the two loop RGEs for SSB masses is given by ||T7 

dt Wn^^""' (167r2)2^™?' ^ ^ 

where t is the natural log of the scale, i = Qj, Uj, Dj, Lj and Ej, and j = 1 — 3 is a 
generation index. Two loop terms are suppressed relative to one loop terms by the square 
of a coupling constant, plus an additional factor of IGtt^ in the denominator. The two loop 
terms 



9 aigja^ + bigla2 + Ciglai, (3) 



where 



ai = \gl{Kml^ + "^L) + Tr[Tnl + 3mi + Sm^, + 2ml + 6m|;]}, (4) 
^2 = gli^H^ + m]j^ + Tr[3mJ + m^]}, and (5) 
^3 = glTr [2m^ + + m^] , (6) 

and the are squared mass matrices in generation space. The numerical coefficients a^, hi 
and Cj are related to the quantum numbers of the scalar fields, but are all positive quantities. 
Thus, incorporation of multi-TeV masses for the first and second generation scalars leads 
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to an overall positive, possibly dominant, contribution to the slope of SSB mass trajectories 
versus energy scale. Although formally a two loop effect, the smallness of the couplings is 
compensated by the much larger values of masses of the first two generations of scalars. In 
running from Mqut to M^^ak, this results in an overall reduction of scalar masses, which 
is felt most strongly by the sub-TeV third generation and Higgs scalar masses, and indeed 
leads to the constraints found in Ref. |T^. For values of SSB masses which fall short of the 
constraints of Ref. , a sort of see-saw effect amongst scalar masses occurs: the higher the 



value of first and second generation scalar masses, the larger will be the two loop suppression 
of third generation and Higgs scalar masses, until the constraint of Ref. |]T3[ takes effect (or 
the lightest SUSY particle (LSP) ceases to be charge or color neutral). 

An example of this effect is shown in Fig. |l], where we plot in frame a) the evolution 
of third generation SSB scalar masses from a common GUT scale value of mo (3) = 1000 
GeV, down to the weak scale. We also take mo(l) = 15 TeV, mi/2 = 1400 GeV, Aq = 0, 
tan /3 = 3 and /i > 0. The dashed curves represent the case for universal scalar masses, with 
mo = 1 TeV for all scalars at Mqut, while the solid lines indicate the GSIMH model. We see 
that at scales close to Mqut, the trajectories of SSB scalar masses differ radically from the 
mSUGRA case, due to the dominant two loop RGE contributions, and that furthermore, 
these contributions overcome the positive one loop contributions from gauge interactions, 
and actually suppress the scalar masses relative to the mSUGRA case with universality. At 
lower energies, the one loop gauge contributions to squark masses again become dominant — 
due to increasing gluino mass and SU{3) gauge coupling — resulting in an upward turn 
of the corresponding mass parameters. However, the final weak scale values of SSB mass 
parameters in this case are suppressed by almost a factor of 2 relative to the model with 
universality. In fact, most of the GSIMH model SSB masses are at sub- TeV levels, as 
opposed to the multi-TeV scale of SSB masses from the mSUGRA model. This means the 
GSIMH model could be in accord with naturalness constraints, even though the mSUGRA 
model is not. A further interesting feature is that, since the two loop RGE contributions 
to m| are larger than those to m|;, the lightest third generation slepton is most likely to 
be dominantly left-handed, instead of n^f/ii-handed, as in the mSUGRA model. In frame 
b), we see that the evolution of SSB Higgs masses is also strongly altered by two loop RGE 
contributions. The absolute values of m% and m% are diminished relative to the mSUGRA 
model case, again resulting in the GSIMH model being more natural. 

The corresponding weak scale sparticle and Higgs mass spectra are shown in Table | for 
the two cases shown in Fig. Two additional cases for a high value of tan /3 = 35 are 
also shown. For the GSIMH model in case 1, the first and second generation scalar masses 
are all ~ 15 TeV, which should be large enough to suppress various flavor changing and 
CP violating processes, with the possible exception oi K — K system. Meanwhile, most 
of the third generation scalar masses are at sub- TeV values, in accord with naturalness 
considerations. The relatively high gluino mass {rrig ~ 3.3 TeV) means that the gluino is 
effectively decoupling as well as first and second generation scalars. It is a general feature of 
viable GSIMH models that not only the gluino, but the other chargino and neutralino masses 
are constrained to be heavy. [| This has important consequences for searches for GSIMH 



We recognize that the electroweak charginos and neutralinos have direct couplings to the Higgs 
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models at colliders. For instance, a Next Linear Collider (NLC) e^e~ machine would need at 
least 1.5 TeV in the CM frame to significantly access the lighter sparticles in this spectrum. 
The last entry in the Table shows the overlap between fi and states, so that the lightest 
stau in the GSIMH model is predominantly left handed. At the CERN LHC pp collider, top 
and bottom squark pair production will be the dominant SUSY process. We have generated 
collider events for all SUSY production processes for case 1 GSIMH model using ISAJET. 
We examined the various multi-jet plus multi-isolated lepton plus missing Et signals using 
the standard cuts and SM backgrounds given in Ref. |jT9|. In none of the signal channels 
examined was case 1 GSIMH model visible above SM background at a ba level, assuming 10 
fb~^ of integrated luminosity. We also looked at possible signals in multijet + events with 
two tagged 6-jets. We required the cuts of Ref. [|0|: -^t > 100 GeV, prib-jets) > 50 GeV, 
and + J2 Erijets) > 1500 GeV (where the sum runs over non-tagged jets). In addition, 
the highest 6-jet was required to have px > 100 GeV. Again, for this case, no signal was 
visible above SM background for 10 fb~^. Alternatively, examination of the "effective mass" 
distribution has been advocated as a means to quickly establish the presence of a SUSY 
signal against SM backgrounds [O]. The effective mass is defined as 



Meff = EtUI) + Et{32) + Erijd,) + ^t(j4)+ ^t, (7) 

where jl refers to the highest Et jet in the event, and so forth. In Fig. |^, we plot the 
effective mass signal for case 1 (open circles) and backgrounds taken from Ref. [ETI], after 



using the cuts of ||2T]]. For the five mSUGRA case studies of Ref. [^, the signal always 
emerges from the background Mgjj distribution at a suitably high value of Me//; for our 
GSIMH case 1, this does not happen. For the spectra of the GSIMH model case 1, clearly a 
more clever dedicated set of cuts will be needed to establish a SUSY signal. Finally, we note 
that the spectra of GSIMH case 1 may well be excluded by constraints from the cosmological 
relic density of neutralinos , since dominant neutralino annihilation likely takes place via 
i-channel stau, sbottom and stop exchange, and these particles are rather heavy. However, 
since m^^ is getting close to nif^ and m^^, co-annihilation effects may be important, and 
could work to reduce the relic density to acceptable levels p3| . 

In case 2, with tan/3 = 35, we take mo(l) = mo(2) = 10 TeV, mo(3) = 900 GeV, 
mi/2 = 1000 GeV and Aq = 0. Again, the third generation scalars are generally suppressed 
to sub- TeV values compared to their mSUGRA counterparts. The top squarks and staus 
are the lightest third generation scalars, with the fi again predominantly left handed. We 
also examined GSIMH model case 2 for visibility at the LHC just as for case 1. Again, we 
found no observable signal using the simple cuts of Ref. [l^ or for the special fe-jets cuts 
listed above. The Me// distribution is shown in Fig. ^ As with case 1, signal is always 
below background levels. The spectra of GSIMH case 2 may well be allowed by constraints 
from the cosmological neutralino relic density, since m^^ ~ and the neutralinos can 
annihilate efficiently through the very wide s-channel Higgs graphs to bb final states. 

To map out the allowable parameter space range for viable GSIMH models, we per- 
formed a scan over parameter space, generating random samples of input parameters over 
the following ranges: 



bosons, and could result in a need for fine-tuning if they are too heavy. 
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mo(l) = mo(2) : 2000 ^ 25,000 GeV, 
mo(3) : ^ 2000 GeV, 
mi/2 : ^ 2000 GeV, 

Aq : -10, 000 GeV ^ +10, 000 GeV, 
tan : 3 55, and 
/i : + or — . 

In Fig. H, each dot represents a viable model. Models are ruled out if i.) REWSB does not 
occur, a.) masses for the physical scalar eigenstates become tachyonic, or in.) the Zi is not 
the LSP. In all frames, the value of mo(l) is plotted on the vertical axis. In a), we plot models 
versus mo(3). It is seen that larger values of mo(l) require larger values of mo(3). In fact, this 
is just the result of Ref. |T^. In frame b), we plot mo(l) versus mi/2- Our main point here 
is that for large enough values of mo(l), there is a minimum value of mi/2 that is required: 
large values of gaugino masses give large negative slope contributions to SSB scalar mass 
trajectories, which can off-set the large positive two loop contributions refered to earlier. In 
general, the lower bound on mi/2 leads to rather large gaugino masses at the weak scale, 
which makes detection of SUSY at colliders more difficult, even though third generation 
scalars can remain at sub-TeV values. In frame c), we plot mo(l) versus the parameter Aq, 
and see that wide ranges of the Aq parameter give rise to viable GSIMH models. Finally, in 
d), we plot mo(l) versus tan/3, and see that no range of tan/? is particularly favoured. 

In Fig. ^, we show sparticle masses and /i contours in the mo (3) vs. mo(l) plane, where 
we take mi/2 = 500 GeV in frames a) and b), and mi/2 = 1000 GeV in frames c) and d). 
We take Aq = and /i > in all frames, and tan /5 = 3 in a) and c) , and tan (3 = 35 
in b) and d). In the dark shaded regions fi^ < so that REWSB fails, while the light 
shaded regions have a non-neutralino LSP; the intermediate shaded region has tachyonic 
scalar masses. In unshaded regions viable models may be possible. Some authors may then 
impose further constraints from fine-tuning, FCNC or other considerations. The constraints 
that we delineate are the minimal constraints any such model should satisfy. We also show 
mass contours of 1 TeV for fi, hi, ii and the fi parameter. The labels are located in the 
regions with mass less than 1 TeV. An exception occurs in frame b), where fi is always less 
than 1 TeV; there we plot a contour of fi = 300 GeV. In frames a) and b), with the smaller 
values of mi/2, the magnitude of mo(l) is generally limited, especially if one requires sub- 
TeV third generation scalar masses, which occur for lower values of mo (3). Moving to higher 
mi/2 values in c) and d), significantly larger values of mo(l) are allowed, which translates to 
greater suppression of FCNC and CP violating processes. Even so, sub- TeV third generation 
scalar masses are still possible in significant parts of the parameter plane. In the large tan P 
plots, we see that mo (3) is restricted from above by REWSB. 

In Fig. we show potentially viable regions in the mo(3) vs. mi/2 plane, for mo(l) = 
mo(2) = 10 TeV, Aq = 0, fi > and a) tan/5 = 3 and b) tan/3 = 35. Points with open 
circles are excluded by REWSB, those with stars by tachyonic masses and finally those with 
dots have a non-neutralino LSP. In addition, various sparticle mass contours are shown. All 
contours show sparticle masses of 1 TeV, except the g contour shows nig = 3 TeV, and the 
Zi contour shows m^^ = 500 GeV. We see that for a fixed value of mo(l) , low values of 
mo (3) give rise to sub-TeV stau masses, while low values of mi/2 give rise to sub- TeV top 
and bottom squark masses. The lightest neutralino is gaugino-like throughout the frame of 
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Fig. ^a, which in general can lead to large values for the relic density of neutralinos, except 
for regions where co-annihilation effects are important (along the boundary of the excluded 
region) . 

At this point, some comparison with the results of Arkani-Hamed and Murayama (AM) 



and Agashe and Graesser (AG) of Ref. |jT3|] seems worthwhile. As an example, taking 



mo(l) = 10 TeV, and mi/2 = 500 GeV, then AM (Fig. 2) require mo(3) > 1500 GeV, 
AG (Fig. 4a) require mo(3) > 1100 GeV and we require, from Fig. |^ for tan/? = 3, 
mo (3) > 2400 GeV. Taking instead mi/2 = 1000 GeV, AM and AG require mo (3) > 300 
GeV, while we require mo (3) > 600 GeV. The work of AM neglects all Yukawa couplings 
in the RGE evolution, and requires only positive SSB scalar squared masses; the work of 
AG is similar, but takes into account the top quark Yukawa coupling. In our work, we 
include the top, bottom and tau Yukawa couplings in the two loop RGBs. Our constraints 
are somewhat different, also, since we require no tachyonic physical masses, a neutralino 
LSP, and REWSB, while the authors of Ref. |]TB|| include some fine-tuning constraints. We 



remark that the exact location of the boundary of the tachyonic excluded region can be 
very sensitive to the manner in which the superparticle mass spectrum is calculated. For 
instance, using ISAJET, the tachyonic excluded region occurs during the first pass of RG 
evolution, when sparticle masses are calculated at tree level. The top quark Yukawa coupling 
can be significantly larger during the first iteration than the second, since loop corrections 
are not yet included. This can cause a greater suppression of top squark squared masses 
than in a fully consistent one loop treatment, for which the bounds on mo (3) would be 
somewhat lower. In a fully consistent one loop calculation of sparticle masses, the excluded 
region should always come from the LSP and REWSB constraints, instead of from tachyonic 
masses. 

Summary and Conclusions: We have investigated supersymmetric models with a scalar 
IMH in place at the GUT scale. Two loop RGEs for SSB masses are crucial for this analysis. 
We map out parameter regions with first and second generation scalar masses in the 5-20 
TeV range but with sub- TeV third generation scalars. These models go a long way towards 
solving the SUSY flavor and CP problems, while remaining on the edge of naturalness. At 
the very least, all models of this type will have to satisfy the minimal constraints that we 
have imposed. Finally, since the GUT scale gaugino mass is bounded from below, charginos, 
neutralinos and gluinos are generally rather heavy, making the resulting sparticle mass 
spectra challenging to discover at planned future collider facilities. 
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TABLES 



TABLE I. Weak scale sparticle masses and parameters (GeV) for two cases of mSUGRA and 
GSIMH models. 



parameter 


mSUGRA 


GSIMH 


mSUGRA 


GSIMH 




case 1 


case 1 


case 2 


case 2 


mo(l) 


1000.0 


15000.0 


900.0 


10000.0 


mo (3) 


1000.0 


1000.0 


900.0 


900.0 




1400.0 


1400.0 


1000.0 


1000.0 


^0 


0.0 


0.0 


0.0 


0.0 


tan P 


3 


3 


35 


35 


nig 


2965.1 


3291.6 


2181.7 


2395.4 




2745.6 


14993.0 


2076.5 


10014.8 


mi 


2639.9 


15009.0 


1997.4 


10024.2 


1351.8 


14983.5 


1110.7 


9991.0 




1122.3 


14991.4 


970.9 


9994.9 


mo. 


1349.9 


14983.4 


1107.8 


9990.7 




2131.2 


658.7 


1616.2 


792.6 




2546.4 


1113.7 


1877.8 


1072.2 




2527.2 


922.3 


1834.7 


940.0 




2638.8 


1316.6 


1907.6 


1147.2 




1121.0 


671.6 


835.4 


690.0 




1351.4 


851.0 


1063.7 


761.2 




1349.4 


667.9 


1055.3 


721.4 




1110.7 


1097.0 


785.3 


745.2 




1110.6 


1097.7 


785.2 


745.8 




600.3 


616.5 


424.6 


432.1 


mh 


112.8 


111.8 


122.4 


126.5 


rriA 


2302.3 


1447.7 


1218.3 


801.9 




2305.2 


1450.7 


1222.5 


807.7 




1721.8 


1200.2 


1070.0 


792.3 




0.02 


0.99 


0.15 


0.71 
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FIG. 1. Evolution of o) third generation SSB masses and b) Higgs SSB masses from Mqut to 
Myjeak versus scale choice Q. Dashed lines indicate mSUGRA model while solid lines indicate the 
GSIMH model. We take mo(3) = 1000 GeV, m.1/2 = 1400 GeV, ^0 = 0, tan/3 = 3 and //, > 0. In 
the GSIMH model, mo(l) = mo(2) = 15,000 GeV, while in mSUGRA, mo(l) = mo(2) = mo(3). 



11 




Meffl(GeV) 



FIG. 2. A plot of the effective mass distribution background (histogram), and signal (open 
circles) for GSIMH model case 1. 
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FIG. 3. A plot of the effective mass distribution background (histogram), and signal (open 
circles) for GSIMH model case 2. 
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FIG. 4. A plot of parameter space where viable GSIMH models arc located. The vertical axis 
is mo(l), while in a) models are plotted versus mo(3), in 6) versus mi/2, in c) versus Aq and in d) 
versus tan /3. 
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FIG. 5. A plot of the ?tio(3) vs. mo(l) plane where viable GSIMH models are located. In a) and 
b), we take mi/2 = 500 GeV, while in c) and d, we take mi/2 = 1000 GeV. In a) and c), tan/3 = 3, 
while in b) and d, tan/3 = 35. We take = and > 0. The contours show where various 
sparticle masses and the /i parameter arc equal to 1000 GeV (except frame b), where /i = 300 
GeV). The darkest region is excluded by lack of REWSB, while the lightest excluded region has a 
non-neutralino LSP. Tachyonic sparticle masses are generated in the intermediate shaded regions. 
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FIG. 6. A plot of parameter space where viable GSIMH models are located in the mo(3) v.s. mi/2 
plane, where we take mo(l) = mo{2) = 10 TeV, Aq = 0, n > and a) tan/? = 3 and 6) tan/3 = 35. 
Points with circles have no REWSB, points with stars have tachyonic masses and dots represent 
points with a charged and/or colored LSP. All contours represent sparticle masses of 1 TeV, except 
nig < 3 TeV below the g contour, and m~ < 500 GeV below the Zi contour. 
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